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Abstract 


Objective: The lifetime difference in adjacent parallel structure components becomes small as the number of components 
belonging to the same parallel structure increases. To infer the system structure, we must clarify the components that 
belong to the same parallel structure. 

Methods: A strengthened change point detection model (SCPDM) for weak mean difference data (WMDD) is established, 
which usually indicates that, as affected by a large variance, the mean difference in two subsignals for one data sequence 
becomes nonsignificant. For repeatedly retrievable WMDD, we performed two enhanced operations that doubled the mean 
difference by using the variance information and analyzed the asymptotic properties of the enhanced data. Then, we 
proposed an SCPDM based on the asymptotic results. 

Results: Finally, we compared the SCPDM with two other main change point detection models and verified that the 


SCPDM is superior to other models using WMDD change point detection by the simulation method. 


Limitations: This paper also have several limitations. First, we only discussed that y,,..., yr are independent with 


normal distribution and single change point. Second, the reason why the relationship between | H-U 3) and O has an 
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important influence on the accuracy of change point detection is not discussed in depth. We only defined the ratio 
boundary of WMDD by experience and simulation. 

Conclusions: Traditional change point detection models may become insensitive or ineffective for WMDD. We gave 
some asymptotic analysis and established a enhanced change point detection model (SCPDM) based on the asymptotic 
results. Compared with the traditional method, SCPDM can effectively detect the change point. 


Keywords: single change point detection; weak mean difference; large variance; enhanced operations; simulation 


1. Introduction 

Weak mean difference data (WMDD) are a type of data in which the information of the mean 
difference is reduced by a large variance. For example, as is shown in Fig 1, when the mean difference 
and the standard deviation , the location of the change point is easily detected. If the standard 
deviation , the accuracy of change point detection may be decreased. However, if the standard 
deviation , the location of the change point is hardly detected by most current models. 

An important example of WMDD comes from the reliability field. Assume that the resistance 
values of all of the components in Fig 2. are equal. Jin et al. [46] pointed out that components which 
are tested in a laboratory environment differ in significant ways from those which have experienced 
operations in fielded systems where the homogeneous components will suffer different degrees of 
damage. As is shown below, the voltages of components belonging to adjacent parallel structures 
become more similar as the number of components belonging to the same parallel structure increases, 
which means that the lifetime difference in components belonging to adjacent parallel structures 
becomes small too. To infer the system structure, we must clarify the components belonging to the 
same parallel structure. As we all know, components in the same parallel structure have homogeneous 
life data. Hence, detecting the small lifetime difference in components is extremely important to 
distinguish whether components belong to the same parallel structure or not and is beneficial to 
establish a topology diagram of the system structure. 

The difficulty of detecting change points for WMDD is how to capture the small differences 
between subsignals. To solve this problem, we performed two enhanced operations which increased 
the mean difference between subsignals by utilizing the information of the variance. In addition, we 
then analyzed the asymptotic property of the enhanced data. Next, we proposed a strengthened change 
point detection model (SCPDM) according to the asymptotic results. Finally, we compared the 
SCPDM with two current main models and verified that the SCPDM has a higher efficiency than 


those of the other models in the WMDD change point detection by the simulation method. 


2. Literature review 

The so-called change point is the location at which a certain amount in the signal sample 
suddenly changes [1]. The change point often represents a qualitative change for our object of focus. 
Historically, Page [2-3] first proposed the study of change points in sample testing. To detect the 
change points in a signal sample, , the process roughly involves the following several steps. First, 
select an associated cost function [4] to measure the homogeneity in each subsignal. Second, 
according to whether the number of change points is fixed, compute a discrete optimal problem to 
obtain the estimations of the locations of the change points. In different change point detection 
models, establishing a suitable cost function for a specific sequence is the first and most important 
step [5]. Many classical change point detection models have been proposed for various kinds of 
signals, mainly including the following three types. For piecewise independent identically distributed 
(i.i.d.) signals, the mean shift model was first established for normal random signal samples which had 
a piecewise constant mean and constant variance [3, 6-9]. Second, certain signals may have their 
means shift along with shifts of their variances. For example, the mean shift and scale shift models 
were established for normal random signal samples with piecewise constant means and variances by 
certain predecessors [10-11]. Except for normal random samples, the rate shift model has been studied 
for Poisson distribution signals with piecewise constant rate parameters [12-13]. The second type of 
change point model is appropriate for a signal with a linear dependency between the variables along 
with changes happening at certain unknown instances, which are also called structural changes [14- 
16]. In this situation, several well-known models were established, such as the autoregressive model 
[1, 17] and multiple regression models [16, 18]. Other commonly used change point detection models 
include kernel change point detection [19-22] and the Mahalanobis-type metric [23]. Kernel change 
point detection can be operated on the high-dimensional mapping of the original signal that is 
implicitly defined by a kernel function. Certain machine learning techniques may be involved in this 
kind of method, such as a support vector machine or clustering [24-25]. In addition, in certain 
clustering methods, the Mahalanobis-type metric is usually used to replace the cost function in the 


mean shift model [26-28]. 


Except for the models mentioned above, several classical models based on algorithms have been 
proposed about inferring change points, mainly including the following four types [10]. The first 
model is based on the likelihood ratio. Csórgó and Horváth [29] established a change point detection 
model under the assumption of a multivariate Gaussian distribution. This model is mainly used in 
analyzing the change points in time series data. The second model is based on the Bayes model. A 
number of researchers have focused on this method. Kander and Zacks [30] aimed at the exponential 
family to establish a change point detection model, while Gardner [31] established a model based on 
the normal distribution. Later, the model was extended to the large sample distribution theory, 
multivariate normal distribution and general linear regression field [9, 32-33]. The third model is 
based on the maximum likelihood. This kind of method has the mature large sample theory. For 
example, Fotopoulos et al. [34-35] established the exact computable expressions, bounds and 
approximations for certain analysis results. The last model is based on samples. This kind of method 
focused on the nonparametric method which has the advantage of being distribution free. 

Among the subdirections of change point detection, signal samples with mean shifts have always 
been a research hot spot. Hawkins et al. [36] used the sample variance without degrees of freedom to 
determine the change point. In [37], the maximum likelihood estimation method was utilized to 
analyze the change point under the premise of verifying the type of population distribution. Later, 
prior knowledge was incorporated in establishing the change point detection model in [38]. In [39], 
the change point location was determined by analyzing the local information near a point, which 
involved complex distribution information usually substituted for by certain approximate results. In 
recent years, certain new methods also have been discussed regarding change point detection. As 
indicated in [40], an optimal algorithm was introduced to determine the location of a change point. In 
[41], an adapted algorithm was established by the polynomial maximization method. In [42], partition 
models were set up for testing the existence of a mean shift and estimating the location of the change 


point. In addition, lasso methods were established and improved by many authors [43-45]. 


3. Methodology (Design/Approach) 


3.1 Two enhanced operations 


Because the standard deviation is far larger than the mean difference, it is unwise and inefficient 
to perform change point detection on WMDD using traditional models. In fact, compared with the 
information of the mean difference, the variance is very remarkable which may supply more 
information. Therefore, we considered utilizing the variance that belongs to a kind of disturbance 


information in analyzing change points by performing the following two enhancements. 


Now, we believe that t° (1^ € {2,..., T3) is the only abrupt location in sequence Jj... Yr . For 
Jis Yr , we first conducted the first operation called enhanced-l at ¢ (t € 也 ;5 了 y) and obtained 


the enhanced-| sequence y, m 
yy = y Y; =min(y, 1, yi) i Ls t- 1l, 


y, =),, y, = max(»,,,),),¡=t+1,...,T. 


Then, we conducted the second operation called enhanced-ll at £ (£ € Q2... Tj ) and obtained 
the enhanced-ll sequence y acids Yo . 


y, 2 y, Y; = max(y; ,,y;) i =1,..,t—1, 


y, = y, Yi = min(y, y), is tels T. 
Intuitively, these kinds of enhancements utilize the variance information directly through taking 


large or one smaller between adjacent samples. 


Now, we provide several symbolic explanations. For a signal sample Jj... Yr, we performed 


the above two operations at location 1 (i € PaT} ), where 44 (f, y ) indicates the sample mean of 


Visio Jis 1, (t, y ) indicates the sample mean of Ypres Yr; 1, (t y”) indicates the sample mean of 


yp Ya; and u;(t,y') indicates the sample mean of Yir Yr; u (t,y") indicates the sample 


mean of Y,»... V, ,. In addition, » and y represent operations on the homogeneous signal sample 


and the signal sample containing a change point, respectively. 


3.2 Asymptotic property of the enhanced sequence 


We considered establishing a certain asymptotic property of the enhanced sequence in the 
following theorem. First, we will explain several symbols. Assume that Y,, with (上 =1,...,7 ), are 
independent random variables, with probability distribution f (v, |0), and@ is a vector-valued 
parameter. Vi...) Vr is a set of samples of Y... Yp. Fora signal sample that has one change location 


at f^ (i € {2,...,T}), 44 indicates the population mean of the first subsignal and 4, indicates the 


population mean of the second subsignal. c^ indicates the constant population variance. Generally, 


WMDD indicate that | Hy “Hal <1. 
O 


Theorem 1. Assume that Yi; is ii.d., i.e., homogeneous with a constant variance o°. The 


above two operations are performed ”, times independently at location 1 (1 € 12,..., Tj ). Then, we 


have the following asymptotic property: 


. 1 n. 
d. im M 


uy (y) 一 ba ey) m lh) 一 ua a?) = | eL 
Proof. Yi». Yr is t.i.d. According to the law of large numbers, we have the following: 


Hi (t, y ) — Has as. [—, 


where 4min is the smaller population mean of the enhanced sequence. 


Likewise, we have the following: 
My; y ) > s a.s. T-t>0 


where Lg is the larger population mean of the enhanced sequence. As a result, we have the 


following: 
Hy SY) a 5) > as LT 
On the other hand: 
us (Ly) Um 4$ to 
uc (Ly) Uu, as. T-t>0 
Consequently: 


Hii (y )- uy (t, y) > Hmax — Hwi QS. t,T -t —> %0 


In addition, f (x,y ) = E 一 y | is a continuous function, so the q.s, converge can be preserved 


under the transformation of this function, so we have the following: 


H min H max 


uy Gy )- uy (6, y") — 


uy (y )- uy (y )2 


Himax ~ min =0 as. tT -t — oo 


Thus, because the process is independently performed M, times, 


D (hy )-3 (69) -名 (y )- ia yn) G =1.....ņ_) can be viewed as independently and 


identically distributed. According to the law of large number, we have the following: 


c i=l 


] (^. | o | 2 . 2 . 
Epa (t, y )- Uy CE- (t, y )7 Hy (t, y joe as. LT -t — on, — 00 


Remark 1. Theorem 1 demonstrates that when Y;.....Y; is homogeneous and the above two 


operations are performed 2. times independently at any location 1 (f € (2,..., T3 ), the value of 


1 n. M pes ， Bs E E Ñ 
(A (t y ) - uy; (£, y | ha, (y ) - Hz GY f) fluctuates at approximately 0. 


c i=l 


Next, we will give the asymptotic property when there is only one change point among the signal 


samples. We will present the asymptotic results at the position of the change point. 


Theorem 2. Assume that the independent V;.... Yr only has one change point £^ with a constant 


variance o? . At £”, the above two enhanced operations are independently performed 2. times. Then, 


we have the following asymptotic property: 


(1) If [A “Hl > 1, we then have: 
O 
is lato... x t ee Le Ro t ues a D DTE TE 
1. EX» Haii (t Y )7 Hai (t y J- Hii (t My )- lxi (t A He. (t Y) 7 i (t y) Hii (t y ) = Usi (t >Y )-o)a, 


The meanings of uy (ty )> dis Vy) , ji (C y) and iis (C y ) are the same as in Theorem 


(2) If LH Hal < 1, we then have: 
O 


1 n. 
P lm 一 
nat ,Tt >o n, ^ 


Proof. (1) First, we prove the first part. The following formula is clear when 4 < LH, 


hana E) 


us (6, y) wef) = 2u, -a =l 


B a) ia a] ju na a) >am) as. tÜ,T-t' o 


| Hy -H| 


Because — — — > 1, an incorrect enhanced operation cannot reverse the direction of the mean 
O 


difference, so we have the following: 


+ uy (y) pu, (€, y") > (u, - u) as. t T-t — 00 


[€ y) 
Likewise, when 44 > H, we have the following two relationships: 


hu a) ha | e a a | n i) as. t',T-t' 0 


" [y )- Ha Co) CF )- Ha Cr a - uj) as. TTi — 
In summary, we have the following: 
十 


uy (y )- us CY) i Ya 6 9) 0 as. C,T -£ >0 


heise] 


iy yt Cy 


Thus, because the process is independently performed n, times, 


neue un) (i=1,...,n,) can be 


+ 


uy (0, y) ua (6, y) 


his a) 1 


Aia ey) 


viewed as being independent and identically distributed. According to the law of large numbers, we 


have the following: 


uy (sy) uy (Cy) HCP wy 0) as. C,T -t',n, >% 


AS el 


ra 
Proof. (2) The following formula is clear when 44 < 45, 


hay Id) - je) Ha (y) > (e - uu) a.s. pem PEO 


| 4 -Hl 


Due to ———— <1, the enhanced-ll operation can reverse the mean difference, so we have the 
O 


following: 


hu € o) is o) - pu ons a) > Ap) 141) as. ¢,T-t +0 


Likewise, when Hı > H2, we have the following two relationships: 


bu a) E- us oa € | de 1) as. tÜ,T-t' — o 


a (t, y") - uy Ek, Cy € a) Sp) as. t,T-t >o 
Thus, 
OA pe oia on] s — as. Ü,T t — o 


Because the process is independently performed n. times, 


- ) can be viewed as being 


soe Tbo 


Ci Cy G=] 


Hy OY )— us Cy) = 
independently and identically distributed. Therefore, we have the following: 


| 


c i-l 


sy) ku C- My (OY) Ha «| > Ju, -u| as. nt T -t > o 


Remark. 2. For WMDD, theorem 2 demonstrates that when YV;....Yr has a change at Ü ( 


te {2,...,7}) and if the above two operations are performed 7。 times independently at /' ( 


1 Ze E ' * ' is , * ý * a * " j or * " 
) the value of 2| Hy © y )- ug (t Y (ty )- My © oy E 


te {2,..., T) nN, al 


reach the maximum value 2|u, —H,|; if the above two operations are performed M, times 


independently at any nonchange location, the value of 


Hii (ty) — Ui eo) will be less than the maximum value 


1 Ne e , » : . " A i 
Al li (t y )= Hoi (t iy ) = 


c i-l 


2/1, = uil . The latter finding is verified in Section 3 in detail. 


3.3 SCPDM 


For WMDD, the different asymptotic properties in theorem 1 and theorem 2 are important 
information for judging whether there is a change point and the location of the point. Consequently, 


we proposed a model to detect the change point in WMDD by sampling repeatedly in this section. 


Assume that there is only one change point whose location is £' (1<t € T) in Yj.-- V7. To 
obtain a better estimate of 三 ， we should establish a contrast function [5] to measure the goodness-of- 


fit of the signal. First, at the location of t (f= Za ), the above two enhanced operations are 


performed 2. times. The following contrast function is established: 


hn 3) la y) | 


‘me 1 n. 
t (y= aema S| 


c i=l 
e . 


My uy 
t 


t 


My 7 My 
t 


t 


According to theorem 2 and theorem 3, the value of is expected to be 


large when f` =t" and small when f is not well-estimated. 


Throughout the establishment process of the SCPDM, how to solve this discrete optimization 


problem becomes clear. We only need to calculate 


1 Ne . | , . i : > a ., x 
-* u y )- us cz? 站 (y )- Ay GY | at ( 


c i=l 


t=T 


hu) is y )- 


RT 1 n. . . e . 
t (y= sema 25 Hy (Ey ) - uu Gy JJ Details can be seen in 


c i=l 


the next section. In the next section, we performed several simulation studies to estimate t under a 


normal distribution with various kinds of parameters. 
4. Results 
In this section, we will perform simulations in two parts. In the first part, we will verify the 


correctness of theorem 2 and remark 2. In the second part, we performed several simulation studies 


to estimate £^ under a normal distribution with different parameters and compared the SCPDM with 


two current main models, verifying that the SCPDM has a higher efficiency than those of the other 


models for WMDD. 


4.1 Verifying the correctness of the theorem 


To verify the correctness of theorem 2, we generated random numbers based on the normal 


distribution, and the parameter settings are shown in the caption of each figure. For the public 


parameters, we set T =1000, 1* = 501, and 7t, 21000, 


| 4, -Hl 


For == > I, as is shown in Fig 3, we carried out two sets of simulation analyses. 
O 


Le “Hal oy 
O 


To verify the correctness of remark 2, when at the location of j ( 


5 


t = 400,...,600 n, 53 


WEDE jeie] is calcúlated and 


the results are shown in Fig 5. 


It can be seen from Fig. 5 that when f is near ¢ =SOl, 


1 n. e. a A e. F " 
| Be (by )— Hi GY | reaches a maximum value, which is 


c i=l 


¿y i y) = 


approximately 2|u, z 1] . The simulation result is consistent with theorem 2. 
5. Discussion 


When LA Hal <1, in order to verify that the SCPDM's estimation of Á is better than those of 


[03 


the other models, we present certain simulated results based on traditional models, including the least 


squares model [47] and Bayes method [48]. We generated random samples V;.... yr based on the 


normal distribution, and the parameter settings are shown in the corresponding figure; we set the 


public parameters, namely, T 21000, ¢ = 500, and 2. = 1000. For the three models with the 


same parameter settings, we repeated the same operation 1000 times, counted the estimation results of 


t^ and regarded the frequency of each Ê ”as the probability of becoming a real change point, which 


reflected the accuracy of each model. The results are shown in Fig 6, Fig 7 and Fig 8. 


By setting different parameters, the detection accuracies of the change point interval, i.e., 


Pr(^€[499,503]), and change point location, i.e., Pr(f” — 501), of the three models were 
compared in Table 1 and Table 2, respectively. 


| 45 -H| 


When ————— >1, ie, the parameters were set to and 


for 1000 
c c -0.5 lau at s 


repeated tests, the Bayes model resulted in Pr(£^ e[499,503]) = 0.87 and Pr(£^ =501)=0.61, 


and the least squares model resulted in Pr(?” e[499,503]) = 0.942 and Pr(f = 501) = 0.561. 
Both methods had a high accuracy for change point detection. 


a Hal eq 
O 


Whe , i.e., the type of data is WMDD and we consider Pr 人 人 € [499,503]) ' 


When the parameters are set to o =| and | Hi - [|= 1 , the accuracy of the SCPDM is 31% higher 


than the accuracy of the least squares model and 41.1% higher than the accuracy of the Bayes model. 


When the parameters are set to o =1 and | H4 -uj- 0.7 , the accuracy of the SCPDM is 54% higher 


than that of the least squares model and 60.7% higher than that of the Bayes model. When the 


parameters are set to o =1 and | Hi - [|= 0.4 | the accuracy of the SCPDM is 62.6% higher than the 


accuracy of the least squares model and 70.1% higher than the accuracy of the Bayes model. When the 


parameters are set to o =1 and | Hı -H IF 0.1 , the accuracy of the SCPDM is 42% higher than that 


of the least squares model and 43% higher than that of the Bayes model. When the parameters are set 


to g =10 and | H4 -uj- 1 , the accuracy of the SCPDM is 31% higher than the accuracy of the least 
squares model and 41.1% higher than the accuracy of the Bayes model. 


| 4, -Hl 
n —— —— 
O 


Whe <1, we consider Pr 人 üy When the parameters are set to and 


gel 
| LH - 4, (+ 1 , the accuracy of the SCPDM is 25% higher than the accuracy of the least squares model 
and 33.6% higher than the accuracy of the Bayes model. When the parameters are set to O =1 and 


| Hi -U= 0.7 , the accuracy of the SCPDM is 27.7% higher than that of the least squares model and 


32.1% higher than that of the Bayes model. When the parameters are set to g = 1 and | 44 -4;|7 0.4 


, the accuracy of the SCPDM is 23.2% higher than the accuracy of the least squares model and 28.8% 


higher than the accuracy of the Bayes model. When the parameters are set to O =1 and 


| 1, -us[- 0.1, the accuracy of the SCPDM is 9.1% higher than that of the least squares model and 


9.5% higher than that of the Bayes model. When the parameters are set to o — 10 and | Hi -uj- 1 ; 


the accuracy of the SCPDM is 9.2% higher than the accuracy of the least squares model and 9.5% 
higher than the accuracy of the Bayes model. 
6. Conclusion 

This paper focused on detecting the change points of WMDD. We conducted a certain 
asymptotic analysis and established an enhanced change point detection model (the SCPDM) based on 


the asymptotic results. According to theorem 2 (2), the enhanced sequence uses the significant 


> 


variance information so that the weak mean difference increases from m TA to 2|lu, — H5 


which makes the change point easier to detect and increases the accuracy of change point detection. In 
addition, for WMDD, the traditional methods may be improved by purely adding the sample capacity 
to the sequence. While under the premise of having the same amount of data, the SCPDM greatly 
increases the efficiency of change point detection by repeatedly detecting sequences with the same 
data structure. In addition, repeated measurements are possible for the life data of components at the 
same location. Hence, compared with traditional methods, the SCPDM can effectively detect change 


points. Although the accuracy of change point detection has been improved, this paper also has several 


limitations. First, we only discussed that JV;..... yr is independent with a normal distribution and a 


single change point. Second, the reason why the relationship between | 44 -L,| and o has an 


important influence on the accuracy of change point detection is not discussed in depth. We defined 
the ratio boundary of WMDD only by experience and simulations. In a future study, we will extend 


the SCPDM to other distribution types and multiple point detection. In addition, for theorem 2, we 


will reprove the theorem by introducing the relationship between | LH, -H " and o. 
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Figure Legends 


Fig 1. Diagram of the accuracy of change point detection influenced by the variance. (A) 0 = 0.1, 
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Fig 7. Probability of each location becoming a change point in the Bayes model. (A) O = 0.5, 
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Fig 8. Probability of each location becoming a change point in the SCPDM. (A) 6-1, 
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Fig.1 Diagram of the accuracy of change point detection influenced by variance. 


1 
—A 
1 4 


la —A —O— 


1A 


m O 


一 《一 
O 
一 人 一 


Voltage difference 2 times 1.5times 1.3times  ... 1 times 
Fig. 2. Varieties in difference between adjacent parallel structural components 


0.03 - 
5 5 
山 Ww 
0.3 
T 0.4 
100 200 300 400 500 0 100 200 300 400 500 
Replication Num ber Replication Number 
(ao =0.1,| 4 -u,ẹ= 0.5 (b)o =1,|A -u}= 5 
Fig. 3. Under 500 repetitions, value of 


tus (Cy) - uy (iy) 


) 


wesie AS 


bi) 


1 
2 


2 EY a — 
0 100 200 300 400 
Replication Number 


(a) o =1,| 4, -Hl 0.1 


100 200 300 400 500 
Replication Number 


(b) o =1,| 1, -4,1=0.7 


100 — 20 300 40 500 
Replication Num ber 


(d) o =10,| 4, -u; 1 


Fig. 4. Under 500 repetitions, value of 


4 : 
0 100 200 300 400 500 
Replication Num ber 
(© o =1,| 4, -4,/=1 
1 n. 
n, j=l 
0.25 


50 100 150 
Location 


(a) o =1,| 4, -4,/=0.1 


200 


My (€, y) 11, (6, y) — y Ey) 


^ 
2:24 
Fi ^ 
P ~ 
Pe ~ 
d b." 
a is 
al , ] 
0 50 100 150 200 
Location 


(5 o =1,| 44 - 4,15 0.7 


2 2.5 
E ^ 2 
P d 要 
= "d ee — 68 
= ite “| = 
es 5 : Loy 
0.51 he 
0 . . 0 : 
0 50 100 150 200 0 50 100 150 200 
Location Location 
(c) o =1,| 4 - 1,21 (d) o =10,| 4, - 1,21 
0.8 
0.85 


1 
Fig. 5. The value of 一 >， 


c i-l 


Probability 


0 200 400 600 


Location 


(a) o 20.5.] 4 -u5]-1 


0.08 


800 1000 


Probability 


Lai 
0 200 400 600 800 


Location 


(d) o =1,| 44 -u= 0.4 


1000 


n 


c 


Probability 


Probability 


450 


500 
Location 


(e) o =1,| 4, -Hl 0.4 


550 


My (by )- hey) - 


location. 
0.4 
03 
0.2 
04 
0 
0 200 — 400 — 600 800 1000 
Location 
(b) o =1,| 4 -u,|=1 
0.2 
0.15 
0.4 
0.05 
0 
200 400 600 800 1000 
Location 


(e) o =1,| 4, -u= 0.7 


Hu (Gy ) - uy (y )| | at 


0.025 


0.01 


Probability 


200 400 600 


Location 


(e) o =10,| 4 - 4,51 


0.02 


800 1000 


0.015 


Probability 


0 200 400 600 


Location 


800 1000 


( o =1,| 4, -1,/50.1 


Fig.6. Probability of the each location becoming a change in Least Square Model 


0.7 = 0.2 0.02 
0.6 
0.5 0.15 0.015 
2 > > 
3 04 Z 5 
E T E 04 E 0.01 
n. o n. 
0.2 0.05 0.005 
E | 
0 0 J 0 
0 200 400 600 800 1000 0 200 400 600 800 1000 0 200 400 600 800 1 
Location Location Location 
(a) EOS fy, po all (b)-o-1,| ui, -u 41 (0.0 10.) y, uu 
n INN an ANO e ^70 jd enn enn inna 
0.06 T T - - 
0.05 ot 0.025 
=| T 
UM > 0.08 002 
Bon $us fus 
Lo 0.04 = 00 
i B B l| Mill 
LLLI | | | ile Tal 
0 200 400 $00 800 1000 0 200 400 600 800 1000 
awi Location Location DM 
dd) of =1,| 4 -H= 0,4 (e)o =1,| 4, SH gp 1,| 4, -4,50 
Fig.7. Probability of the each location becomi g a change in Bayes Model 
0.1 
0 1 1 1 , 0 
200 400 600 800 1000 0 200 400 600 800 1000 
Location Location 
(c) o =1,| 4 “HF 1 (d) & =10,| 4, -u}= 1 
0.35 ; 
0.3 
0.25 
> 
= 0,2 
a 
oO 
RARE 
n. 
0.1 
0.05 
0 L L 1 L 
0 200 400 $00 800 1000 
Location 
(e) & =1,| 4, -41,/=0.4 
Fig. Probability of the each location becoming a change in SCPDM 
格式 要 求 


; 图 注 和 表 注 : TNR, 8 pt. 


is 


图 字 、 表 字 ; 图 题 、 表 


KR: 单 栏 最 大 : 185 mm x 120mm; WERK: 宽 不 超过 60 mm. 


N 


201904.00096v1 


chinaXiv 


图 清晰 度 : 300 dpi 或 以 上 。 若 图 中 有 多 个 图 表 ， 用 A、B、C 等 分 别 表 示 。 


注 明 图 〈 横 纵 坐标 ) 表 的 度量 衡 单位 及 其 说 明 。 


按照 正文 中 引用 的 顺序 升序 排列 ， 每 个 图 单独 占据 一 页 ， 


AN: 


公式 放 在 正文 中 ， 需 用 公式 编辑 器 进行 编辑 。 


不 同 图 之 间 上 月 


Tables 


Table.1 Estimate probability of a change in certain intervals of three models. 


Change point location 


Model — | 44, -4 = 484-488 489-493 494-498 499-503 504-508 509-513 514-518 
0.5 0 0.001 0.028 0.942 0.029 0 0 
Least 1 1 0.010 0.031 0.113 0.681 0.113 0.021 0.013 
— 1 10 0.006 0.008 0.010 0.019 0.013 0.007 0.012 
quars 0.4 1 0.054 0.075 0113 0215 0.107 0.055 0.043 
Neth 0.7 1 0.039 0.048 0.166 0.417 0.148 0.048 0.034 
0.1 1 0.009 0.014 0.015 0.010 0.015 0.011 0.009 
1 0.5 0 0 0.07 0.87 0.06 0 0 
1 10 0.01 0.01 0.01 0.01 0 0 0.01 
Bayes 1 1 0.01 0.02 0.11 0.58 0.13 0.07 0.03 
0.7 1 0.06 0.06 0.18 0.35 0.18 0.07 0.01 
Method 0.1 1 0.01 0.03 0.01 0 0.02 0 0 
0.4 1 0.02 0.04 0.08 0.14 0.10 0.06 0.02 
1 1 0 0 0.001 0.991 0.008 0 0 
SCPDM 1 10 0.019 0.065 0.171 0.430 0.201 0.073 0.023 
0.7 1 0 0 0.010 0.957 0.033 0 0 
Method 0.1 1 0.019 0.065 0.171 0.430 0.201 0.073 0.023 
0.4 1 0 0 0.040 0.841 0.118 0.001 0 


Table. 2. Estimate probability of ( = t“ of three models 


MD SD Least Squared Method Bayes Method SCPAM 
1 0.5 0.651 0.61 --- 
1 1 0.286 0.2 0.536 
1 10 0.003 0 0.095 
0.4 1 0.066 0.01 0.298 
0.7 1 0.164 0.12 0.441 
0.1 1 0.004 0 0.095 
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